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GROWTH DUE TO BUOYANCY OF WEAK HOMOGENEOUS TURBULENCE WITH SHEAR 

by Robert G. Deissler 
Lewis Research Center 

SUMMARY 

An analysis is made to determine how the energy in a weak homogeneous turbulence 
with buoyancy and shear changes with time. It is shown that, although the turbulence 
ultimately decays when buoyancy is absent, the presence of buoyancy counteracts the 
decay, and the turbulent energy grows at large times. 


INTRODUCTION 

Studies of weak homogeneous (grid generated) turbulence with uniform mean velocity 
gradient are given in references 1 and 2. In the results obtained there, the turbulent 
energy always decayed with time. The energy produced by the mean velocity gradient 
was at all times less than that dissipated. In reference 3, where the initial condition 
was modified to give a finite initial turbulent energy, the energy sometimes increased 
for a while, but it still ultimately decayed. This behavior was attributed to the fact that, 
while the total energy was increasing, energy was being drained out of the component of 
the turbulence in the direction of the velocity gradient by the pressure-velocity correla- 
tions. Since there was no turbulence production in that component, it quickly decayed, 
with the result that the turbulent shear stress, and consequently the turbulence produc- 
tion in all the components, ultimately decreased. 

Thus, the key to obtaining a nondecaying turbulent shear flow appears to lie in 
keeping the energy from being drained out of the turbulence component in the direction 
of the mean velocity gradient. In strong turbulence, the distribution of energy among 
the directional components is evidently accomplished by the pressure-velocity correla- 
tions, but in weak turbulence, those correlations generally tend to make the turbulence 
more anisotropic (ref. 1). If, however, we superimpose on the shear flow destabilizing 
buoyancy forces in the direction of the mean velocity gradient, it may be possible to 



obtain a nondecaying solution, even for weak turbulence. Those buoyancy forces should 
tend to prevent the turbulence component in the direction of the mean velocity gradient 
from decaying. 

Some of the effects of combined buoyancy and shear on weak homogeneous turbulence 
were investigated in reference 4. However, the question of whether or not the turbulence 
decays was not considered. It is with that aspect of the problem that we will be con- 
cerned here. 


ANALYSIS 


The basic equations for the present study can be written as (ref. 4) 
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where the subscripts (except e) can take on the values 1, 2, or 3, and a repeated sub- 
script in a term signifies a summation. The quantity u^ is a fluctuating velocity com- 
ponent, t is the time, U- is a mean velocity component, p is the density, p is the in- 
stantaneous pressure, p g is the equilibrium pressure, x. is a space coordinate, v is 
the kinematic viscosity, j3 = -(l/p)(ap/3T) , g. is a component of the body force per 

.H 1 

unit mass, t is the temperature fluctuation, T is the mean temperature, a is the 
thermal diffusivity, and overbars designate averaged values. (Symbols are defined in 
the appendix. ) 

In obtaining equations (1) to (3), the instantaneous velocities and temperatures in 
the incompressible Navier-Stokes and energy equations are broken into mean and fluc- 
tuating components. For the last term in equation (1) (buoyancy term), the density is 
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assumed to depend effectively only on temperature and is not far removed from its equi- 
librium value (value it would have for no heat transfer or turbulence). The equation for 
the pressure (eq. (2)) is obtained by taking the divergence of the Navier-Stokes equation 
and applying continuity. 

Equations involving correlations between fluctuating quantities at two points P and 
P' can be constructed from equations (1) to (3). Such equations have been obtained for 
homogeneous turbulence with uniform velocity and temperature gradients in reference 4. 
| They can be converted to spectral form by taking their Fourier transforms. We con- 

sider the case where the velocity and temperature gradients are in the x^-direction 
^ (vertical) and the body force (gravity) is in the -x^-direction. Let g = -g^, a = dUj/dXg, 

and b = dT/dx^, and assume that the turbulence is weak enough to neglect terms con- 
taining triple correlations. The correlation equations then become in spectral form 
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The quantity k is the wave number vector and d/T = dx^ d/Cg d/Cg. The magnitude of ~k 
has the dimension 1/length and can be considered to be the reciprocal of a wavelength or 
eddy size. The unprimed and primed quantities in the barred products refer, respec- 
tively, to values at points P and P’ separated by the vector r. The quantity 6^ is 
the Kronecker delta. 

Equations (4), (5), and (6) give contributions of various processes to the rate of 
change of spectral components of uTuT, tuT, and t 2 , respectively. The second term in 
each equation is a transfer term which transfers activity into or out of a spectral com- 
ponent by the stretching or compressing of turbulent vortex filaments by the mean veloc- 
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ity gradient, as discussed in reference 1. The terms with k in the denominator are 
spectral components of pressure-velocity or pressure-temperature correlations and 
transfer activity between directional components (ref. 1). The terms proportional to 
/3g and 6^ (or 6jg) are buoyancy terms which augment or diminish the activity in a 
spectral component by buoyant action. The last terms in the equations are dissipation 
terms, which dissipate activity by viscous or by conduction effects. The remaining 
terms in the equations produce activity by velocity or temperature gradient effects. 

For solving equations (4) to (6), the turbulence is assumed to be initially isotropic 
at t = tp. That condition is satisfied by 
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where Jq is a constant that depends on initial conditions and Kq is an initial wave 
number that is characteristic of the turbulence. Equation (10) differs from the initial 
condition in references 1, 2, and 4 (but not from that in 3) by the exponential factor, 
which was set equal to 1 in those references (kq = °°). For the initial conditions on 6 
and y. (at t = t^), it is assumed that 

6 0 ' ( r i) 0 = 0 < U > 

That is, the turbulence producer (grid) is assumed to be unheated, so that the tempera- 
ture fluctuations are produced by the interactions of the mean temperature gradient with 
the turbulence. 

A method of solving the preceding set of partial differential equations (eqs. (4) 
to (6)) is described in reference 4. Those equations are converted to ordinary differ- 
ential equations, and spherical coordinates are introduced by using the transformations 
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Kj = k cos <p sin 9, k 2 ~ K s ^ n s ^ n and K % = K cos The resulting set of equations 
is integrated numerically by machine computation. Directionally integrated spectrum 
functions can then be obtained from 
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These spectrum functions can be integrated over all wave numbers to give the following 
single-point correlations: 
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Computed correlations and spectra for the case where the buoyancy forces are desta- 
bilizing (negative vertical temperature gradients) will be considered in the next section. 
The results given there are for a gas with a Prandtl number of 0. 7. 


RESULTS AND DISCUSSION 

The effect of destabilizing buoyancy forces on weak homogeneous shear-flow turbu- 
lence in a gas is illustrated in figure 1. The superscript (a) on uTuT^, t^ a \ and 
indicates that those parameters have been made dimensionless by using quantities re- 
lated to the shear (in contrast to those related to the buoyancy, which will be used later). 
Curves are shown for two values of Richardson number and of the initial wave number 
parameter. 

The curves indicate that for a Richardson number of 0 (no buoyancy effects) all com- 
ponents of the turbulent energy decrease with time. The turbulent shear stress -u 1 u„ 

fal 1 6 

also decreases with time, except near the initial time. (At t v ; = 0 the turbulence is 
isotropic and UjU^ is 0.) 
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(b) Dimensionless initial wave number, 




ul/2* 

(dU 1 /dx 3 ) 1/2 


Figure 1. - Effect of destabilizing buoyancy on 
variation with time of weak turbulence in uniform 
uniform shear flow. Prandtl number, 0.7. 


The decay of the components of turbulent energy for no buoyancy effects evidently 
occurs mainly because there is no production term in the equation for u^u^ (component 
in the direction of the mean velocity gradient). This can be seen by letting i = j = 3 in 
equation (4), in which case the production term (first term on the right side) drops out. 
In addition, the pressure-velocity correlation terms in equation (4) tend to drain energy 
out of the u|-component when the turbulence is weak, as discussed in reference 1. As 
a result the u|-component decays rapidly compared with the other components, which 
have energy fed into them by the mean velocity gradient or by the pressure-velocity 
correlations (see fig. 1). When u| decays, the shear component u^u^ must also de- 
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cay. There is then no mechanism for maintaining the turbulence since that maintenance 
apparently takes place as a result of work done on the turbulent shear stress by the 
velocity gradient. (See first term on right side of eq. (4).) All of the turbulence com- 
ponents must then decay. 

By contrast, for Ri = -1 (buoyancy forces destabilizing), all components of the tur- 
bulent energy decay for a while and then begin to increase without limit as time becomes 
large. This increase evidently occurs because the vertical buoyancy forces excite the 
u|- component of the turbulence and replenish the energy being drained out of it. 

It might seem surprising that all components of the turbulence continue to increase 
with time rather than level off. There are no boundaries on the flow considered here 
however, so that the effective Reynolds number and Rayleigh number of the mean flow 
are infinite. As the scale or mixing length of the turbulence continues to grow, the 
eddies encounter larger and larger velocity and temperature differences, so that the ef- 
fective driving forces on the turbulence continue to grow. 

Comparison of figures 1(a) and (b) shows, as expected, that for p) = «> (all wave 
numbers present), the components of the initial energy are infinite, whereas for = 1 
they have a finite value. The turbulent shear stress -UgUj starts at zero on both plots 

since the turbulent shear stress for isotropic turbulence is zero. For the case of 
( a \ 

k o - °°, however, the value of -u^u^ jumps to infinity in an infinitely short time and 
then decreases. For = 1, (- u^u^) first increases steadily and then either decreases 
(Ri = 0) or continues to increase (Ri = -1). 

To give an idea of the distribution of the turbulent energy with wave number, energy 
spectra (spectra of urn. ) are plotted in figure 2 for Kq & ) = °° and Ri = 0 and -1. For 
P) = 0, i is proportional to k ^ (eq. (10)). As time increases, the spectra move to 
the smaller wave number regions; that is, the scale of the turbulence grows indefinitely 
large with time, since the fluid is unbounded. 

Thus far we have been considering the effect of buoyancy on a shear-flow turbulence. 
Next we want to consider the related problem of the effect of imposing a mean shear on 
turbulence that is buoyancy-controlled. For doing this it is convenient to use the param- 
eters uf^), p), Kq 3 ), T*0°) t and Tup), which have been made dimensionless by using 
quantities related to the buoyancy (see figs. 3 and 4). The parameters used in figure 1 
were, on the other hand, nondimensionalized by using quantities related to the shear. 

The effects of shear on buoyancy- controlled turbulence are illustrated in figure 3 

and 4, where u?^ 3 ), t2^), and Tup) are plotted against p) for several values of 
1 /wv 

Richardson number and k ^ . For the case of no shear (Ri = -°°) the results were ob- 
tained from the integrated equations in reference 5. All components of the turbulent 
energy, as well as the temperature fluctuations and the temperature-velocity corre- 
lations, increase as p) becomes large. This occurs even when shear is absent and the 
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(b) Richardson number, Ri = -1. 


Figure 2. - Variation with ti me of turbulent 
energy spectra (spectra of UjUj). Prandtl num- 
ber, 0.7. 
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Dimensionless time, = (3^^g^ 2 (dT/dX3) 1 ^(t - 1 q) 
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Figure 3. - Effect of uniform shear on variation 
with time of weak turbulence in flow with 
destabilizing buoyancy. Prandtl number, 0.7. 
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Dimensionless time, t (b) = (J 1/2 g 1/2 (dT/dx 3 ) 1/2 (t - t Q ) 
(b) Dimensionless initial wave number, 

(b)K (b). .! 

0 MT/dx/VV /4 

Figure 4. - Effect of uniform shear on varia- 
tion with time of temperature fluctuations 
and temperature-velocity correlations in 
flow with destabilizing buoyancy. Prandtl 
number, 0.7. 
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turbulence is completely controlled by the destabilizing buoyancy forces (Ri = 

Although all turbulent energy components can increase with time when shear is absent, 
the component in the direction of the buoyancy forces is, in that case, at least an order 
of magnitude greater than the other components. On the other hand, when both buoyancy 
and shear are present, all components can be of the same order of magnitude. 

Loeffler has recently considered the effect of a gradient in electric charge and an 
applied electric field on homogeneous turbulence (ref. 6). That problem is analogous 
to the present one for the case when no shear is present and the Prandtl number 
is infinite. It was found that the turbulent energy increases without limit as time in- 
creases, when the electric field is in the direction of increasing charge density. For 


large times the turbulent energy was proportional to (exp t^)/^^) . 

As the shear increases (as Ri goes from to -2), the turbulent activity in 
general increases, at least at the earlier times. The shear does not seem to affect 

or 7Ug^ at the smaller times when = 1, however. At larger times, although 

Uj and u| increase with increasing shear, u|, r^, and tu^ all decrease with increas- 
ing shear. These decreases appear to be related to the fact that at large times the 
presence of the shear causes energy to be drained out of the u| component (as dis- 
cussed earlier), and thus out of tu^ and t 2 (see eqs. (5) and (6)). 


CONCLUDING REMARKS 

Although weak homogeneous turbulence with a uniform shear ultimately decays with 
time, the presence of destabilizing buoyancy forces in the direction of the mean velocity 
gradient prevents that decay. In that case the buoyancy forces replenish the energy 
being drained out of the component of the turbulence in the direction of the mean velocity 
gradient, and the turbulent energy increases without limit as time increases. Appar- 
ently the energy can increase without limit because the effective Reynolds and Rayleigh 
numbers are infinite in an unbounded fluid. As the scale or mixing length of the turbu- 
lence continues to grow, the eddies encounter larger and larger velocity and tempera- 
ture differences, so that the effective driving forces on the turbulence continue to grow. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 6, 1970, 

129-01. 
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APPENDIX -SYMBOLS 


a 

vertical velocity gradient, 

n 

defined by eq. (8) 


dUj/dXg 

A 

defined by eq. (12) 

b 

vertical temperature gradient, 

6 

defined by eq. (9) 

g 

dT/dx 3 

vertical body force/unit mass 

5 ii 

Kronecker delta 

in -x^-direction (gravitational 

0,<P 

angular coordinates 


force/unit mass), -g 3 

K i 

wave number component 


body force component /unit mass 

K 0 

characteristic initial wave 

g 3 

vertical body force/unit mass 

number (see eq. (10)) 

J 0 

constant that depends on initial 

V 

kinematic viscosity 

conditions 

P 

density 

Pr 

Prandtl number, v/a. 

T 

temperature fluctuation 

P 

pressure 

TU i 

temperature-velocity correlation 


P g dT/dx 3 

<Pii 

defined by eq. (7) 

Ri 

Richardson number, 




(dUj/dXs) 2 

\L.. 

defined by eq. (12) 

r, r i 

vector between points P and 

Subscripts: 


P' 

1 

in flow direction 

T 

mean temperature 

3 

in vertical direction, which is 

t 

time 


direction of mean velocity 

t 0 

U i 

u i 

x i 

initial time 

mean velocity component 
fluctuating velocity component 
space coordinate 

gradient and buoyancy force 
Superscripts: 

' at point P’ 

- average 

(a) parameter nondimens ionalized 

a 

thermal diffusivity 


by quantities related to shear 

P 

expansion coefficient 

(b) 

parameter nondimensionalized by 

r. 

7 

defined by eq. (12) 


quantities related to buoyancy 
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